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^^ Inelastic light scattering spectra of several members of the RFe3(B03)4 family reveal a cascade 

.^^ ' of phase transitions as a function of temperature, starting with a structural, weakly first order, 

phase transition followed by two magnetic phase transitions. Those consist of the ordering of the 

r^ ' Fe-spin sublattice revealed by all the compound, and a subsequent spin-reorientational transition for 

Cn ' GdFe3(B03)4. The Raman data evidence a strong coupling between the lattice and magnetic degrees 

of freedom in these borates. The Fe-sublattice ordering leads to a strong suppression of the low 

energy magnetic scattering, and a multiple peaked two-magnon scattering continuum is observed. 

D ' Evidence for short-range correlations is found in the 'paramagnetic' phase by the observation of a 

X I , broad magnetic continuum in the Raman data, which persists up to surprisingly high temperatures. 



PACS numbers: 75.40.Gb (Dynamic properties dynamic susceptibility, spin waves, spin diffusion, dynamic 
scaling, etc.) 63.22.-|-m (Phonons or vibrational states in low-dimensional structures and nanoscale materi- 
als) 75.30.-m (Intrinsic properties of magnetically ordered materials) 64.70.p (Solidsolid transitions) 78.30.-j 
(Infrared and Raman spectra) 



I. INTRODUCTION 



The family of RM3(B03)4 crystals (R=Y or a rare earth, M=A1, Ga, Sc, Cr, Fe) has attracted considerable attention 
during the last years. Different combinations of R and M lead to a large variety of physical properties that, together 
^ I with their excellent physical characteristics and chemical stability, make these crystals extremely interesting both 
lO ' from application and fundamental points of view. The lack of inversion symmetry (they crystallize in a trigonal space 
'""' , group) has stimulated different applications in the field of optical and optoelectronic devices. Crystals of YAI3 (803)4 
22 ' ^iid GdAl3 (603)4 doped with Nd have been widely studied and used in optical devices, for self-frequency doubling 
^r \ and self-frequency summing laser aAi^ . Concentrated NdAl3 (603)4 crystals are efficient media for mini-lasers^. Rare- 
vQ ■ earth borates with magnetic ions (M — Cr, Fe) are much less studied. Recently, some iron borates were reported 
f^ ■ to possess multiferroic features, demonstrating the coexistence of magnetic and ferroelectric order parameters^; they 
"^ . may therefore be considered for optoelectronic applications. These magnetic borates are also expected to exhibit 
d ' interesting magnetic properties because of the presence of two different kinds of magnetic ions (3d and 4f elements), 
G , and in particular because their structure features isolated helicoidal chains of magnetic 3d atoms. 
I The crystal structure of RM3 (603)4 borates belongs to the structural type of the natural mineral huntite 

' O ] CaMg3(C03)4 that crystallizes in the space group i?32 of the trigonal system^*Si2i&. The primitive unit cell con- 
tains one formula unit. Three kinds of coordination polyhedra are present, trigonal prisms for ROe, octahedra for 
MOg, and two types of planar triangular BO3 groups: equilateral BIO3 and isosceles B2O3 (the numbering of inequiv- 
alent ions position is the same as in Ref 1^) . The MOg octahedra share edges forming mutually independent helicoidal 
chains which run parallel to the c-axis (see Fig^. The ROg prisms are isolated polyhedra, each of them connects 
three helicoidal MOg chains, while there are no direct R-O-R links between different ROg prisms. This structure can 
';J! ■ also be seen as a succession of alternating BO3 and R-l-M layers perpendicular to the c axis. 
5y , Recently, a systematic study of rare-earth iron borates RFcs (603)4 (R — Y, La-Nd, Sm-Ho) was undertaken for 

polycrystalline samples^. The measurements of magnetization, specific heat, and ^^Fe Mossbauer spectra revealed an 
antiferromagnetic ordering at low temperatures. The magnetic ordering temperature increases with decreasing R'^"'" 
ionic radius, from 22 K for LaFes (603)4 to 40 K for TbFe3 (603)4. In addition to this, X-ray diffraction, specific heat 

[width=Omm]figl.eps 

FIG. 1: (Color online) The structure of RM3(B03)4 incorporates helicoidal chains of M''"'" ions. The picture shows the structure 
of iron borates (M=Fe) viewed from different angles a) perpendicular to the C3 axis and b) parallel to it. The possible exchange 
paths are shown. 



measurements, and differential thermal analysis have indicated a structural phase transition for RFea (603)4 com- 
pounds with R = Eu-Ho, Y. Its temperature increases with decreasing R'^^ ionic radius, from 88 K for EuFca (603)4 
to 445 K for YFe3 (603)4. Magnetic ordering of the La- and Nd-iron borates has been confirmed by spectroscopic and 
magnetic measurements on non oriented single crystal o^i^" . Studies on oriented single crystals have only been carried 
out for GdFe3(B03)4, using a variety of methods^iSiiiii^iiiii. Recent specific heat, Raman scattering, and Nd^+ probe 
absorption measurements revealed a cascade of three spontaneous phase transitions in GdFe3(B03);^. A structural, 
weakly first-order, phase transition into a less symmetric low-temperature phase was found at Ts=156 K (this in 
contrast to the powder result T5=174 K^). At lower temperature, two magnetic transitions have been observed. A 
second-order phase transition occurs at 37 K, resulting in an antiferromagnetic ordering of the Fe spin sublattice. This 
transition is accompanied by a polarization of the Gd spin subsystem. The second magnetic transition is observed at 
9 K, and has been identified as a first-order Fe spin-reorientation transition^^. A detailed picture of magnetic phase 
transitions in GdFe3(B03)4 as a function of temperature and magnetic field has recently been suggested in Ref.^'^ on 
the basis of antiferromagnetic resonance studies. 

Quite recently, the low-temperature structural phase of GdFes (603)4 has been identified as a P3i21 trigonal 
structure*. To our knowledge, no information exists on the low-temperature structure of other rare-earth iron borates. 
Such information can be obtained from the comparison of low-temperature Raman spectra of GdFe3 (603)4 and 
other compounds from the RFe3 (603)4 family. Raman studies of RFe3 (603)4 were performed for R=La and Nd at 
room temperature-^. Preliminary temperature dependent results on the R=Gd compound have also been reported 
recentlyii^. 

The present work investigates two different aspects of the rare-earth iron borates. First, inelastic light scattering 
studies of the vibrational and structural properties of different RFe3 (603)4 are presented. Amongst others they reveal 
the weak first-order phase transition from the high-temperature R32 structure to the low-temperature P3i21 one. 
The second part of this paper focuses on the magnetic properties. Apart from observing changes in phonon frequencies 
at the antiferromagnetic ordering transition arising from magneto-elastic couplings, the spectra also exhibit relatively 
strong two-magnon scattering features. The magnetically ordered phases exhibits well defined magnon excitations 
resulting in a multiple peaked scattering continuum at finite energy. In the paramagnetic phase a strong low energy 
continuum is found, which results from the persistence of short range correlations in the helicoidal Fe spin chains up 
to quite high temperatures. 

II. EXPERIMENTAL DETAILS 

Different samples of RFe3(603)4, with R — Gd, Nd, Tb, Er, and Y were grown using a 6i2Mo30i2 - based flux, as 
described in RefiiSiii. Unhke the 6i203 based fluxesi, in 6i2Mo30i2 - based flux 6i203 is strongly bonded to M0O3 
which excludes a partial substitution of bismuth for a rare earth during crystallization. Spontaneous nucleation from 
the flux resulted in small single crystals ( 1x1x1 mm). We used these crystals as seeds to grow large ( 10x5x5 mm) 
single crystals of Gd, Nd, and Tb iron borates. As for ErFe3(603)4 and YFe3(603)4, only small single crystals of 
these compounds were studied in this work. All the crystals were green in color. The samples were oriented either by 
X-ray diffraction technique or by using their morphology and optical polarization methods. 

The Raman measurements were performed in a backscattering configuration, using a three-grating micro- Raman 
spectrometer (T64000 Jobin Yvon) equipped with a liquid nitrogen cooled charged coupled device (CCD) detector. 
The frequency resolution was better that 1 cm^^ for the frequency region considered. The samples were placed in 
a optical microscope cryostat. The temperature was varied from 2.7 to 500 K, with a stability of ±0.02 K. The 
scattering was excited by the second harmonic light of a Nd:YV04 laser (532 nm), focused down to '--^50 fiin^ with 
the power density on the sample kept below 0.1 mW//Ltm^. The polarization was controlled both on the incoming and 
outgoing beams giving access to all relevant polarization combinations. 

Most of the experiments were performed with light incident perpendicular to the C3-axis of the crystal, fc _L c (as 
a rule, precise orientation of k in the (a6)-plane was not known), however, when the shape of the samples allowed, 
also the k\\c configuration has been employed. Spectra are marked using the Porto notation!^. For k J- c, only 
the polarization of incident and scattered light are indicated. The conventions z||C3, x||C2 are used to describe the 
geometry of the Raman experiments. We were able to obtain only partially polarized Raman spectra of ErFe3 (603)4 
andYFe3(603)4- 

Preliminary dielectric measurements were performed for RFe3 (603)4 on non oriented crystals of typical size 
1x1x0.2 mm^. The surfaces were polished and covered with Ag paste and the electrical contacts were made us- 
ing Pt wires connected to the surface by additional Ag paste. The samples were mounted on a home-made insert for 
a Quantum Design PPMS system. The capacitance was measured using a Andeen-Hagerling AII2500A capacitance 
bridge operating at a fixed frequency of IkHz. 
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FIG. 2: Polarized Raman spectra of GdFe3(B03)4 at room temperature. The bars at the top indicate the frequency regions 
for the internal vibrations of BO3 groups. 
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FIG. 3: (a, b, c) Polarized Raman spectra for RFe3(B03)4, R=Gd, Nd, Tb, at room temperature. (d,e) Partially polarized 
Raman spectra for (d) ErFe3(B03)4 at 390 K and (e) YFe3(B03)4 at 420 K. 

III. EXPERIMENTAL RESULTS 

The room-temperature (RT) spectra of GdFcs (603)4 for all studied experimental configurations are presented in 
FigEl The regions where the internal BO3 vibrations are expected are indicated by the horizontal bars in the top of 
the panel. The spectra show a strong dependence on the polarization of the incoming and scattered light, but also 
show a distinct difference between (xx+yy+xy) polarized spectra for k parallel to c and for k perpendicular to c. The 
spectra for different rare earth iron borates are, as expected, qualitatively quite similar, as is shown in Fig|31 

The phonon lines narrow progressively with lowering the temperature and several new modes appear abruptly in 
GdFe3(B03)4, TbFe3(B03)4, ErFe3(B03)4, and ErFe3(B03)4, at rs=155 K, 198 K, 340 K, and 350 K, respectively. 
Figures El and and Tab. IIIII give a detailed comparative picture of the phonon modes observed above and below Tg 
in RFc3(B03)4. The intensities of the new modes show a hysteresis as a function of temperature (see FigE)). The 
lowest-frequency intense new mode appears abruptly below the phase transition, but then exhibits a soft-mode-like 
behavior with further lowering of the temperature; i.e. for GdFe3 (663)4 and TbFe3 (603)4 the low frequency mode 
shifts from 26 cm~^ at Ts to respectively 56 cm~^ and 57 cm~^ at 3 K, and for ErFe3(B03)4 and YFe3(B03)4 it 
goes respectively from 25 to 72 cm~^ and from 27 to 76 cm~^. At the same time its linewidth decreases from 12 to 
3 cm~^ (see FigsQand 4b). The spectra in (xz-l-yz) polarization of all the compounds demonstrate a low-frequency 
scattering continuum that gradually changes its shape with lowering the temperature and exhibits the opening of a 
"gap" at Tjv (see Fig|Hl. For NdFe3(B03)4 at the lowest temperatures a well-defined narrow peak at about 10 cm^^ 
is observed within the gap. This mode shifts to zero frequency and strongly broadens upon approaching Tjv from 
below (FiglSfc)). For the Gd compound a similar mode is observed below T/v at about 18 cm~^. 

IV. GROUP-THEORETICAL ANALYSIS 

A. High-temperature structure R32 {D3) 

The primitive cell of the high-temperature structure i?32 of RFe3 (603)4 contains 20 atoms which results in 57 
vibrational normal modes. Knowing the local symmetry of all the atomic positions, one can perform the factor-group 
analysis to find the symmetries of these modeaiSiiS. R and Bl atoms reside in highly symmetric D3 positions^iS. and 
generate A2®E modes each. Fe, 62, 01, and 02 atoms are in C2 symmetry positions, which results in Ai©2j42®3-E' 
modes for each of them. 03 atoms occupy general Ci positions and give 3Ai©3A2©6£' modes. Summing all these 
modes and subtracting the A2 ® E acoustic modes one gets the following optical vibrational modes of the crystal, in 
accordance with the results of Refii^. 

T„br = 7Ai(xx,yy,zz) ® 12^2 (^//z) 

©19^(^//x,£;//y,xz,yz,xy) (1) 

Notations in parentheses refer to the allowed components of electric dipole moment (infrared (IR) activity) and the 
polarizability tensor (Raman activity). Doubly degenerated E modes are polar and both IR and Raman active. The 
light propagating along the c-axis {k\\c) probes TO modes, while the kJ-c configuration gives pure LO modesi^. 
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FIG. 4: Raman spectra of GdFe3(B03)4 above and below 7^ = 156 K, namely, at 161 K (the lowest traces in each panel), 150 
K (middle traces), and 30 K (upper traces), (a) (zx-|-zy) polarization; (b) (zz) polarization. New modes that appear below Ts 
are indicated by arrows. Asterisks mark the lines from another polarization. 
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FIG. 5: Comparison of several low-temperature new modes in different RFe3(B03)4. Spectra in (xx+xy)-polarization for 
R=Gd and Tb and partially polarized spectra for R=Y and Er. 
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FIG. 6: The lowest-frequency Ai-mode appearing in low-temperature P3i21 structure demonstrates hysteresis in a plot of 
its intensity vs temperature. The data for Gd(Tb)-compound are presented by open (black) triangles. Inset shows this mode 
during cooling down for temperatures 154.91, 155.03, 155.08, 155.11, 155.16, 155.20, 155.45 K in Gd-compound. The lower 
temperature axis refer to R=Gd, the upper one - to Tb. 

The strongest bonds in the structure of RFca (603)4 borates are B-0 bonds within the planar BO3 triangles. 
Considering them as molecular units that occupy D3 (for BIO3) and C2 (for B2O3) positions, and adding the modes 
of translational character coming from R and Fe ions (respectively in D3 and C2 position) , we obtain the following 
optical modes: 

r^'Jf, = 2Ai © 5^2 © 7E; T^^^^, ^ Ai ® 3^2 © 4E, (2) 

where the librational ones (right equation) come from the BO3 units. The total number of external optical vibrations 
is, thus, 

Tl^^l; ^ 3Ai ® 8A2 ® UE (3) 

Subtracting ^ from Q), one gets the internal vibrational modes related to the vibrations of BO3 groups: 

r;'|* = 4^1 ® AA2 ® 8E (4) 

The correlation scheme of FigElhelps to get a deeper insight into the origin of the vibrational modes of RFe3 (603)4. 
It follows from this scheme that two of the seven Ai modes originate from the fully symmetric vibration of the BIO3 
and B2O3 groups, two correspond to the B2O3 2^3 and v^ vibrations, and two to the rotational and translational 
modes of the B2O3 group. The remaining Ai mode results from displacements of the iron atom. Table HI indicates 
the origin of crystal modes in the frequency ranges of internal vibrations of the BO3 molecules. 

B. Low-temperature structure P3i21(D|) 

The primitive cell of the low temperature structure contains three formula units® (60 atoms) which gives 177 
vibrational modes. Unique positions for the R and Bl atoms remain but their symmetry lowers from D3 in the i?32 
structure to C2 in the P3i21 one. On the contrary, instead of single threefold C2 symmetry positions for Fe, B2, 01, 
and 02 there are now two different positions, namely, a threefold C2 and a sixfold Ci positions for each of them. 03 
oxygen atoms occupy now three different general (Ci ) positions^. The factor-group analysis of the low-temperature 
structure P3i21 gives the following optical vibrational modes: 

r^,b(P3i21) =27^1 ©32^2© 59^; (5) 

Carrying out the same procedure as in the previous subsection one finds what modes correspond to the external 
(translational and librational) motions of the R and Fe atoms and of the BO3 entities and to the internal vibrations 
of the BO3 "molecules". The results are summarized in Table Hll 

V. DISCUSSION 

A. Phonons of RFe3(B03)4 in the R32 and P3i21 phases 

First of all, we consider the polarized room-temperature Raman spectra of RFe3 (603)4, R=Nd, Gd, and Tb 
(Figl^l). Table UTTl summarizes the observed Raman modes of these iron borates at RT. It has been constructed taking 
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FIG. 7: (Color online) The lowest-frequency mode appears at the phase transition at a finite frequency (characteristic of a 
first-order transition), but exhibits a high-frequency shift with further lowering the temperature (characteristic of a second- 
order transition). The inset shows the derivative of the peak frequency in the vicinity of the magnetic ordering temperature, 
evidencing the coupling of vibrational and magnetic degrees of freedom (see text for discussion). 

into account the results of the group-theoretical analysis for the i?32 structure and the frequency ranges for normal 
modes of the regular planar free ion (BOs)'^^ (they are indicated in Fig|21 see also Tables III and ITl|l . The spectra of 
ErFe3(B03)4 and YFe3(B03)4 were only partially polarized. Comparing the 390 K spectra of ErFe3(B03)4 and the 
420 K spectra of YFe3(B03)4 with the RT spectra of other iron borates we complete Table UTTl for these compounds 
as well (at room temperature, ErFe3(B03)4 and YFe3(B03)4 are already in the low-temperature phase). 

The low-frequency part of the Raman spectrum is dominated by external modes arising from translational motions 
of the R and Fe atoms and from librational and translational modes of BO3 groups. The number and polarization 
properties of the modes observed below 500 cm~^ agree with the group-theoretical predictions. All the modes observed 
in the frequency ranges of the v^, 1/2, and ui vibrations of BO3 correspond well to those derived using the correlation 
analysis and summarized in Fig. |5| In contrast, there is one extra mode in the region of the v^ vibration. Possibly, 
the strong doublet at about 1220 cm~^ in the (zx-l-zy)-polarization arises due to the Fermi resonance^ between the 
1^3 vibrations and an overtone of the 1/4 vibration. 

The highest vibrational frequency (originating from the v^ vibration of the BO3 group) observed in ErFe3(B03)4 
differs markedly from the highest frequencies of both TO {k \\ c) and LO (fc _L c) phonons (1280 and 1414 cm-1, 
respectively, in GdFe3(B03)4). This is because in our particular non oriented ErFe3(B03)4 sample the direction of k 
was arbitrary, while the frequencies of polar E modes depend on the direction of k. 

Summarizing the RT results, it is concluded that (i) the number and symmetries of the Raman modes observed for 
the RFe3(B03)4 in the high-temperature phase are fully consistent with the predictions of group-theoretical analysis 
based on the R32 space group; (ii) the modes above ~550 cm~^ originate from the internal vibrations of the BO3 
groups, while those below ~550 cm~^ can be considered as external modes. 

In the low-temperature P3i21 phase, 20 Ai (12 external and 8 internal) and 40 E (24 external an 16 internal) new 
modes should appear, according to the results of the group-theoretical analysis (see Tablenj. These new modes arise 
due to two reasons. First, the symmetries of the local positions for some of the atoms reduce. The intensity of these 
new modes should be proportional to the deviation from the former symmetric position. Second, the primitive cell 
contains now not one but three formula units, leading to an additional Davydov (factor-group) splitting proportional 
to the strength of interaction between equivalent atoms inside the new primitive cell. Both these effects are, as 
a rule, small. That is why the number of the observed new modes (lOAi © 18£' in GdFe3(B03)4) is lower than 
the number predicted by group-theoretical analysis. The same new modes appear below Tg in all the compounds 
RFe3(B03)4, with R=Gd, Tb, Er, and Y (see Fig. 5 and the bottom of Table III), which strongly suggests the same 
P3i21 low-temperature structure in all of them. 

B. Weak first-order structural phase transition 

The most characteristic Raman feature that announces the phase transition from the high-temperature i?32 struc- 
ture to the low-temperature P3i21 one is the sudden appearance of a strong new low-frequency mode (see Figs01ini 
andEJ in the Raman spectra. As the intensities of new modes are proportional to the squares of atomic displacements, 
/ ~ (5^, the strongest new mode is believed to be associated with the biggest displacements. A detailed analysis, ac- 
cording to Refi^, of the structural changes shows that the biggest displacements are those associated with the BO3 
"molecules". In particular, BO3 triangles, perpendicular to the C3 axis in the i?32 structure, tilt by ~ 7° in the 
P3i21 phase; the Bl atoms shift by ^^0.03 A from the centers of regular triangles. The shifts of boron ions relative to 
neighboring oxygen ions create local dipole moments; their triangular arrangement corresponds to an antiferroelectric 
ordering below the temperature of the structural phase transition. This ordering manifests itself via a strong dielectric 
anomaly at Ts observed in our preliminary dielectric measurements (Fig. I10|) . 

The structural changes considered above give rise to many new Raman active vibrational modes connected with 
the BO3 groups (see Table (nj, in particular, to 4Ai additional librational modes. The energy value of the intense 
excitation measured is in the typical range of the molecular librations. At 3 K, it is 56, 57, 72, and 76 cm^^ for 
Gd(157), Tb(159), Er(167), and Y(89) compounds, respectively. As the values for the Er and Y compounds are 
very close, notwithstanding a big difference in atomic masses of Er and Y, we conclude that the rare earth does not 
take part in this lowest-frequency vibration. A difference between the values for the compounds with relatively big 
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FIG. 8: (Color online) Raman spectra of different RFe3(B03)4 collected in (xz+yz) polarization. The upper picture (a) 
shows the opening of a gap in the magnon Raman scattering at Tjv = 37K for GdFe3(B03)4 (the intensity scale of the left 
bottom corner of the picture is optimized to evidence the low-frequency mode, and the dotted lines are guides for evidencing 
the scattering features). The central picture (b,c) show the magnetic scattering around T^ for YFe3(B03)4 and TbFe3(B03)4. 
The lower one (d,e) shows the two-magnon spectra in the full temperature range (2.7K< T <300K) for NdFe3(B03)4 on the 
left and for GdFe3(B03)4 on the right. For the Gd sample a phonon line is evident in proximity of the magnon feature. LFP 
indicates the low-frequency peak, see the text for the discussion. Above Tn the paramagnon scattering survives up to at least 
nearly ten times the ordering temperature. Different temperature measurements in the panels b,c,d,e are depicted with an 
arbitrary offset for clarity. 

[width=0mm]fig9.eps 
FIG. 9: Correlation scheme for vibrational modes of RFe3(B03)4, space group R32. 

(Gd, Tb) and small (Er, Y) ions comes from the difference in interatomic distances and, hence, force constants. The 
ensemble of these observations strongly suggests that the phase transition in RFea (603)4 is correlated to a rotational 
mode of the BO3 groups. 

The abrupt appearance of new phonon modes and the hysteresis observed for their intensities when approaching T5 
from different sides indicate the first-order character of the phase transition. The energy of the low frequency mode 
strongly increases upon lowering the temperature below Ts (Fig[7|). At the lower temperature (3K) the frequency 
of the mode has more than doubled. Such an anomalously large shift is typical for soft modes that announce a 
second-order structural phase transition. Apparently, we deal with a so-called "weak first-order" phase transition. It 
is worth noting that the lowering of the symmetry from i?32 (Uj) to P3i21 (Z3|) is indeed also compatible with a 
second-order phase transition^^. Thus, the observed first-order character of the phase transition does not depend on 
symmetry changes but rather can arise from the third order therm in the Landau expansion of the crystal free energy 
h, allowed by the crystal symmetry, or from the negative sign of the fourth order coefficient u: 

G = Go + rQ^ + hQ'^ + uQ* + ..., (6) 

where Q is the order parameter—. 

Figure [TTl shows the temperatures of the structural phase transition in RFca (603)4 plotted versus ionic radii of 
R"^^, as determined in Ref.^ for powder samples and in the present study for single crystals. A mismatch between our 
data and those of Ref. ^ could arise from a rather poor quality of the Er and Y iron borate single crystals. A further 
study is necessary to clarify this question. 

Finally we note that NdFca (603)4 preserves the i?32 structure down to the lowest temperatures, while RFcs (603)4 
with R'^''" smaller than Tb'^^ have the low-temperature P3i21 structure already at room temperature. 

C. Magnetoelastic coupling 

The inset of Fig[7| shows the derivative of the peak frequency of the strongest new Ai mode as a function of tem- 
perature in the vicinity of the magnetic ordering temperature T/v . The coupling of vibrational and magnetic degrees 
of freedom is evident. The discontinuity in the phonon frequency observed at T = T^v is ascribed to magneto-elastic 
coupling. The magnetic ordering causes spontaneous magnetostriction, that is, atomic displacements. The latter, evi- 
dently, infiuence vibrational frequencies. Recently, Ref."* reported on the study of magnetoelectric and magnetoelastic 
properties of GdFe3(B03)4. A strict correlation between both has been shown experimentally and discussed theoreti- 
cally. The authors of Ref.^ assumed that, below T/v, the crystal symmetry lowers, due to magnetoelastic coupling with 
spontaneous magnetic moments lying in the ab-plane. This breaks the symmetry of the antiferroelectric arrangement 
of the electric dipole moments in the P3i21 structure and leads to the appearance of an electric polarization which 
could be responsible for the weak growth of the dielectric constant in the temperature region Tr < T < Tj^, (see 
also Fig. I10|l . The proposed lowering of crystal symmetry below T/v, however, should result in the appearance of new 
vibrational modes. No indications for additional modes below T/v have been found. This does not necessarily mean 
that the model by Zvezdin et al. is in error. It might very well be that the structural changes are simply too small 
to be detected by Raman scattering. 

As Reference^ suggests, the spin-reorientation first-order phase transition at Tr into the antiferromagnetic config- 
uration of spins parallel to the c-axis recovers the P3i21 structure (that does not allow the presence of the electric 
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FIG. 10: Dielectric measurements of GdFe3(B03)4. The structural phase transition is made evident by a step in the capacitance 
at Ts- No evidence is found of the Fe ordering phase transition at Tjv, while the capacitance shows a discontinuity at the 
spin-reorientation phase transition at Tr (displayed in more detail in the Inset). 
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FIG. 11: (Color online) Temperatures of the structural phase transitions for different RFe3(B03)4 as a function of the ionic 
radius. Our data on single crystals are shown by circles, those of Refj- on powder samples are presented by squares. 

dipolc moment) and, therefore, is expected to be accompanied by a jump of the dielectric constant. As is clear from 
Fig. 1101 the dielectric function indeed exhibits such a jump, thereby strongly supporting the suggested scenario. 

D. Magnetic scattering 

As discussed in the introduction, due to the structural properties of the RFcs (603)4 family, Fe atoms (S=5/2) 
are arranged in helicoidal chains. As a result, the Fe-Fe distance inside the chain is substantially smaller than the 
distance between different chains. The main exchange interaction is therefore expected to be between Fe ions on the 
same chain {J//))J±), suggesting that low dimensionality could play a crucial role in the magnetic properties of the 
family. This, together with the presence of two coupled magnetic sublattices makes the compounds with magnetic 
R'^+ ions good candidates for exotic magnetic ground states. 

A complex magnetic behavior has indeed been found and widely discussed in GdFe3(B03)4. The first magnetization 
measurements on oriented single crystals were interpreted under the assumption that both Fe and Gd subsystems 
order antiferromagnetically at T/v=40 Kii. In this model, below 10 K, the total magnetic moment of the three Fe 
sublattices (S = 5/2) oriented at polar angles ~ 60 ° to the c axis is compensated by the Gd moment (S = 7/2) 
oriented along the c axis. Sharp singularities in magnetization in the magnetic field H||c observed below 10 K were 
attributed to spin-flop transitions in the Fe subsystem. It was assumed further that in the range of temperatures 
between 10 and 40 K the magnetic moments of every Fe sublattice fall into the ab plane, perpendicular to the c axis, 
preserving the 120 ° azimuthal orientation. 

This model was later reconsidered in Ref.^^. It has been suggested, on the basis of the analysis of the Nd- probe 
spectra and of the speciflc heat data, that the Gd subsystem does not undergo a magnetic phase transition at T/y but 
only gets polarized by the Fe subsystem. The following model has been put forward. At T^—Sl K, the iron magnetic 
moments order antiferromagnetically in the direction perpendicular to the c-axis and polarize the Gd subsystem. Below 
Tfi=9 K, the magnetic moments order antiferromagnetically along the c-axis. The temperature Tr=9 K corresponds 
to a first-order spin-reorientation magnetic phase transition^^. Simultaneously and independently, the same conclusion 
has been drawn in Ref.^? where magnetic structures and magnetic phase transitions in GdFea (603)4 were studied 
using antiferromagnetic resonance experiments. This study revealed also a detailed picture of the magnetic structure 
and anisotropy for different magnetic phases in GdFca (603)4 as a function of both temperature and magnetic field. At 
T/v=38 K, the Fe magnetic subsystem orders into a two-sublattice collinear easy-plane antiferromagnet and polarizes 
the Gd spins, which also form a two-sublattice antiferromagnetic subsystem. The anisotropy constant of the Gd 
subsystem has an opposite sign to that of the Fe subsystem. The Gd contribution to the total anisotropy energy 
grows with lowering the temperature, in conjunction with the growing polarization of the Gd subsystem, and becomes 
appreciable below ~20 K. At T^^IO K the total energy of anisotropy changes its sign which results in the spontaneous 
spin reorientation transition. 

The long-range order of the spins on the iron subsystem manifests itself in the inelastic light scattering experiments 
with two main features. The first one, as reported in FiglSJ is the arising of a low-frequency peak (LFP) at ~10 cm~^ 
in NdFe3(B03)4 and at ~18 cm~^ in GdFe3(B03)4. These peaks soften and broaden upon approaching T/v. The 
energy of 10 cm^^ in NdFcs (603)4 has been identified, by absorption spectroscopy of Nd^+ crystal-field levels, as the 
exchange splitting energy of the Nd"^^ ground Kramers doublet at 5 K arising from the interaction with the ordered 
Fe sublattice^". The temperature dependences of the frequency and the linewidth of the Raman LFP are the same 
as those of the ground-state splitting and, respectively, the ground-level width found from optical spectroscopjsiS. 
Therefore the 10 cm"-'^ Raman mode observed in NdFe3 (603)4 is assigned to spin flip scattering on a single Nd'^+ 
moment in the effective field created by the Fe sublattice. Most likely, the 18 cm~^ excitation observed in GdFe3(603)4 
has the same origin. The absence of LFP in Raman scattering of YFe3 (603)4 with nonmagnetic Y'^+ (see Fig. ^ is in 
favor of such interpretation. The crystal field of RFe3 (603)4 splits the ground-state multiplet of an ion R'^^ with odd 
number of electrons into r4 and T^q Kramers doublets within the D3 site symmetry. The analysis of Raman selection 
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FIG. 12: (Color online) First moment < lu > oi the magnetic excitation vs temperature calculated from 5 to 85 cm~^ for 
GdFe3(B03)4 (circles) and NdFe3(B03)4 (squares). The inset shows the integrated intensity of the magnetic signal at different 
temperature, compared to the magnetic heat capacity of YFe3(B03)4 (solid line)-. 

rules for the ground-state spin flip scattering shows that in the case of the r4 ground state all the components of 
the Raman tensor are allowed while in the case of the T^q ground state only the antisymmetric xy-yx component is 
Raman active. Thus, a relatively strong (xz) Raman LFP of NdFca (603)4 points to the r4 symmetry of the Nd'^"'" 
ground state. 

The second major feature of the iron spin ordering is the arising of a broad structured scattering band around 
60 cm~^ ascribed to two-magnon Raman scattering involving the creation of a pair of magnons with wave vectors k 
and -k. Below the Neel temperature this is a characteristic feature for all the different compound investigated of the 
RFe3(B03)4 family (R=Y, Er, Tb, Gd and Nd). Figure |H1 shows the temperature evolution of this broad signature 
of magnetic scattering. The spectra reported in Fig|Sl are found in the (xz+yz) polarization only. No evidence of 
magnetic scattering is found when the k vector of light is parallel to the c-axis (z(xx-l-yy-l-xy)z). 

As discussed in the previous section, NdFe3 (603)4 does not undergo a structural phase transition and the crystal 
space group remains the high temperature one (i?32), while the symmetry of all the other compounds investigated 
(for R=Gd, Tb, Y, and Er) is reduced to the space group P3i21. It is evident from FiglHlthat the magnetic scattering 
spectra of all iron borates are quite similar. So, the presence of inequivalent Fe chains^ in the low-temperature 
structure of RFe3(B03)4, R = Gd, Tb, Y, and Er does not strongly affect the magnetic excitation spectra, and the 
spectra can be analyzed taking into account only the four different magnetic ions present in the i?32 unit cell (three 
Fe'^"'" ions and one R'^+ ion). Still, a rather complicated two-magnon spectrum is expected, strongly depending on the 
anisotropy parameters. At the lowest temperature reached (T=2.7 K) the two-magnon spectra indeed show a complex 
structure exhibiting at least three main peaks (FiglSJ. As the most efficient mechanism of two-magnon scattering in 
antiferromagnets is usually the exchange-scattering mechanism, and the strongest exchange interaction is between the 
iron atoms along the helicoidal chains, these peaks presumably arise from three magnon branches representing spin 
excitations on the iron chains. 

The broad magnetic scattering feature of Gd, Tb, Er and Y compounds is approximately centered around ~70 cm~^, 
while the one of Nd is centered at lower frequency (~60 cm"-'^). This suggests that the bigger ionic radius of Nd'^+, 
causes the interatomic Fe-Fe distances to be larger, and therefore the exchange interaction to be smaller, this would 
be consistent with the scaling of T/v with the rare earth's ionic radii^. Unfortunately, due to the presence of phonon 
modes at the same frequency of the magnetic excitation, it is not possible to determine the frequency scaling of the 
magnetic scattering among the other compounds (Gd, Tb, Er and Y) to confirm it. 

At the temperature of the spin-reorientational transition (Tr), no drastic changes are observed in the magnon 
spectra. This observation confirms that the observed scattering is mainly due to magnetic excitations on the Fe'^+ 
chains, and sheds light on the proposed exotic magnetic ordering of Fe spins for the GdFe3 (603)4 compound. The 
configuration with 120° angle between nearest neighbor Fe'^^ spins for Tn > T > Tfj which then reorient at Tr 
into an easy axis antiferromagnet along the c-axis as proposed inii is not consistent with our observation. Such a 
drastic change in the spin configuration would induce a major change in the two-magnon dispersion and inelastic 
scattering spectra. The picture of a reorientational phase transition at Tr, as proposed iniS*i^, survives assuming 
an easy-axis/easy-plane anisotropy perpendicular to the c-axis in the temperature region Ta < T < Tn- Indeed, 
considering a simple Heisenberg-type Hamiltonian, where the interaction depends only on the nearest neighbor scalar 
product, an easy axis anisotropy perpendicular or parallel to the c-axis would produce similar magnetic excitation 
spectra and therefore similar two-magnon spectra. 

A striking feature of the magnetic scattering, as shown in Fig. |S1 is the persistence of the magnetic scattering in an 
anomalous shape up to extremely high temperature above T^- In the three dimensional antiferromagnet NiF2 this 
"paramagnon scattering" has been observed well into the paramagnetic phase up to 4T/\^. The first moment of the 
frequency < w > as a function of temperature in NiF2 and MnF2 exhibited a fast decrease to lower frequency (~ 10% 
of the "0 temperature value") at Tn and a subsequent smooth decrease to zero for T > Tm, in excellent qualitative 
agreement with the theory prediction for the moments^i. Also in the case of RFe3 (603)4 the two-magnon spectrum 
softens continuously its frequency approaching the Neel temperature (FiglH)). Yet, surprisingly, the paramagnon 
scattering survives at all measured temperatures above T^ {i.e. up to about 10 times T/v). As reported in Fig. 1121 
the first moment of the excitation decreases at T/v less than expected and approaches zero only very slowly above T/v. 
Considering NdFcs (603)4, the minimum at T/v and subsequent increases with temperature up to 2.5 T/v of the first 
moment is possibly due to the presence of luminescence lines overlapping with features of the magnetic scattering. 
This is likely also the origin of the temperature shift to high frequency (lower absolute energy) of the excitation 
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FIG. 13: (Color online) Low frequency Raman scattering above the Neel temperature for GdFe3(B03)4. The different temper- 
ature are depicted with an offset for clarity. The paramagnon scattering is present up to almost ATn ■ The fit is obtained with 
an fenomenological curve /((jj) = jS^ + L, where G is a Gaussian distribution centered at frequency, and L is a Lorenzian 
distribution to take into account the presence of the vibrational mode (see text for discussion). 

observed at around SOcm^^. 

From a theoretical point of view the presence of quasi-elastic scattering vi^ell into the paramagnetic phase of three- 
dimensional antiferromagnets has been ascribed to two possible origins, namely, spin diffusion^^ and spin density 
fluctuations or energy difFusion^^. In the scenario of spin diffusion, absent in a perfect one-dimensional system, the 
four-spin correlation function (describing the Raman response) is well approximated by a Gaussian shape centered 
at zero energy. In contrast, for the energy diffusion scenario one expects a lorentzian lineshape in the absence of 
strong spin-lattice interactions. The integrated intensity Im in the energy diffusion scenario is expected to follow the 
magnetic contribution to the heat capacity Cm as Im oc CmT'^'^^^. In the presence of strong spin-lattice coupling this 
later relation is expected to remain valid, even though the low energy scattering is expected to broaden significantly 
due to the rapid relaxation of the magnetic excitation energy into the lattice. 

The shape of the low-energy scattering feature in the paramagnetic phase is different from any previously reported. 
Neither a Lorenzian nor a Gaussian distribution can fit our data. This anomalous shape of the paramagnon scattering 
continuum is a common feature of all the compounds investigated; Fig^ldepicts the observed spectra of GdFes (603)4 
as an example. The observation that the spectra can not be modelled with either a Gaussian nor a Lorentzian 
response might evidence the presence of important spin-lattice interactions, in line with the observations discussed in 
the previous section. A comparison of the integrated intensity for GdFea (603)4 and NdFe3 (603)4 to the magnetic 
heat capacity determined by Hinatsu et al^ (see inset fig. 12) for YFe3 (603)4 does not show a fair agreement in the 
'paramagnetic' phase, even though there is a good correlation below the phase transition. Apparently energy diffusion, 
combined with spin-lattice coupling does not explain the observed scattering. On the other hand, as shown in Figure 
El the spectra for GdFe3 (603)4 can be fitted nicely at all temperature with a broad phenomenological function 
consisting of a renormalized Gaussian distribution I{lo) = . ^^l . -t- L, where G{u!) is a pure Gaussian, and L(ijo) is a 
lorentzian to account for the phonon line observed in the spectra. This equation is reminescent of the response function 
of interacting magnons in low dimensional quantum spin systems, although it would be unclear why this would still 
hold for a. S — 5/2 system at high temperatures. At present one can not draw a definite conclusion, but the results 
do strongly suggest that short range order correlations in the spin system survives up to quite high temperatures. 
Most likely an interplay between the low dimensionality and strong spin-lattice interactions are responsible for the 
observed unique behavior of the magnetic scattering up to several times T/v, and deserves further studies. 

VI. SUMMARY 

In summary, structural, magnetic, and magneto-elastic properties of various members of the RFe3 (603)4 family 
(R = Gd, Tb, Nd, Er, and Y) have been studied using primarily inelastic light scattering. The compounds show a 
weak first-order structural phase transition which results in the activation of a strong mode of 6O3 librational nature, 
and its subsequent temperature evolution. The scattering spectra observed in the low-temperature phase are fully 
consistent with the low-temperature structure reported earlier. A detailed analysis of the vibrational spectra at the 
temperature of the magnetic ordering transition showed a strong magneto-elastic coupling for all the compounds. The 
analysis of the two-magnon Raman scattering in the magnetic phases showed qualitatively similar scattering spectra 
for the five compounds investigated, evidencing that the structural differences between them do not strongly affect 
the magnetism and the magnetic excitation spectra. Finally an unprecedented an intriguing paramagnon scattering 
continuum up to quite high temperatures. The origin of this scattering does not seem to be the usual energy diffusion 
observed in a variety of other magnetic system, and seems to indicate the presence of short-range order spin-spin 
correlations arising from the low-dimensional nature of these compounds. 
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TABLE I: Frequency range and assignment of the vibrations of the BOa-group. The second column hsts the assignments for 
the free ions. The third column shows the expected internal modes and their symmetries in the high temperature phase of the 
rare-earth iron borates (the superscripts 1 and 2 refer to the different borate groups in the rare-earth iron borates). 
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TABLE II: The origin of vibrations for two structural modifications of RFe3(B03)4 and the new modes that should appear 
below the temperature of the structural phase transition Ts . 

species T > Tsi R32 (Dl) T < T^: P3i21 (Dj) New modes 

'R 
Fe 

external BO3 
translational 
librational 
internal BO3 
Total 7Ai © 13^2 © 20s 27 Ai © 33^2 © 60S 20Ai © 40E © 2OA2 
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®E 




Ai 


©2A2 


(BSE 


Ai 


©3^2 


ffi4S 


Ar 


©3^2 


ffi4S 
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iffi4^ 


2©8£; 



Au 


5 2^2 ( 


B3E 
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©5A2 


®9E 


5Ai 


©7^2 


®12E 


5Ai 


®7yl2 


®12E 


12A 


.1 © I2A2 © 24E 



Aii 
3Ai 


© 2B © ^2 

© 6£ © 3^2 


4Ai 
4Ai 
8A1 


© 8S © 4^2 
© 8S © 4^2 
© WE © 8A2 
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TABLE III: Observed vibrational modes for the high-temperature phase of several rare-earth iron borates. 
External modes {3Ai ® HE) 



NdFe3(B03)4 
Ai Eto 



180 
298 
473 



89 



260 
272 
312 
354 



475 



Elo 

93 

159 

193 

232 

266 

281 

332 

356 

384 

439 

488 



GdFe3(B03)4 
Ai Eto 

180 84 

307 160 

475 195 

270 
273 
315 
352 
391 
443 



Elo 

93 

160 

198 

232 

270 

287 

330 

357 

391 

443 

488 



TbFe3(B03)4 

Ai E 



ErFe3(B03)4 
Ai E 



180 
309 
476 



93 

159 

198 

230 

270 

287 

330 

355 

392 

444 

480 



181 
312 

472 



84 

160 

199 

225 

272 

325 
350 
395 

442 
475 



YFe3(B03)4 
Ai E 
180 107 
312 158 
468 200 

268 

327 
354 
390 
441 

470 



Internal modes (4Ai ® 8S) 



NdFe3(B0 

Ai 

636 

950 

990 

1220 



3M 



E 

579 

625 

669 

734 

968 

1195 

1218 

1244 

1260 

1413 



GdFe3(B03)4 

Ai Eto 

638 

957 631 

990 670 

1230 



1198 
1280 



Elo 

580 

633 

676 

735 

968 

1212 
1229 

1250 
1414 



TbFC3(B03)4 

Ai E 



ErFe3(B03)4 
Ai E 



635 
957 
988 
1220 



579 
632 
676 
734 
966 

1214 
1230 

1254 
1411 



632 

960 
988 
1230 



579 
632 
675 
730 
960 

1220 
1233 

1342 



YFe3(B03)z 
Ai E 
631 576 



984 



Modes appearing in the low T phase: 20^i (12 cxt. and 8 int.) ffi AQE (24 ext. an 16 int. 



NdFe3(B03 



GdFe3(B03)4 
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53 

144 

203 
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244 

263 



368 
385 
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724 
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114 

167 

206 

254 
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281 

305 

310 

337 

374 

395 
472 
596 
654 
669 
955 
968 



TbFe3(B03)4 

Ai E 

54 



ErFe3(B03 



672 

732 

959 

1219 

1230 

1408 



YFe3(B03)4 



232 
247 
265 



368 



679 
723 
953 



169 
207 
256 
278 
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306 
311 
337 
375 

398 
473 
597 
654 
670 
954 
966 



72 
149 
202 
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E 
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208 
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373 



470 



Ai 


E 
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370 
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